Momentum dependent light scattering in insulating cuprates 
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We investigate the problem of inelastic x-ray scattering in the spin— i Heisenberg model on the 
square lattice. We first derive a momentum dependent scattering operator for the A\ g and B\ g polar- 
ization geometries. On the basis of a spin-wave analysis, including magnon-magnon interactions and 
exact-diagonalizations, we determine the qualitative shape of the spectra. We argue that our results 
may be relevant to help interpret inelastic x-ray scattering experiments in the antiferromagnetic 
parent cuprates. 



The advances made in 3 rd generation light sources have 
recently provided new insights into the study of electron 
dynamics in strongly correlated systems via resonant in- 
clastic x-ray scattering (RIXS). Detailed information has 
already been obtained on Mott gap excitations and or- 
bital transitions as a function of doping in the cuprate 
familiesi 2 -! 3 -^^. However, due to limitations in resolu- 
tion, ~ 300 meV near the elastic line, excitations at low 
energies remain hidden. On the other hand, low energy 
excitations at small photon momentum transfers, e.g., 
phonons, magnons, and electron-hole excitations in met- 
als and superconductors, have been well characterized via 
Raman spectroscopy^. Yet, achieving a detailed under- 
standing of the momentum and polarization dependence 
of low energy excitations in strongly correlated matter 
would greatly clarify the interplay between various de- 
grees of freedom, such as antifcrromagnctism, charge den- 
sity wave order, and superconductivity, often present in 
system with many competing interactions^. 

In the energy range below a few hundred meV lies one 
of the most prominent features observed via Raman scat- 
tering in Heisenberg antiferromagnets - the two-magnon 
feature at energy ~ 2.7J, with J the nearest-neighbor 
magnetic exchange. While the peak frequency of the 
broad two-magnon peak in Raman scattering is well un- 
derstood, the asymmetry of the lineshape as well as 
the polarization dependence remains unexplained, even 
though it has been lavished with attention^. Prelimi- 
nary RIXS on the Cu K -edge in LaCuO^ and M-edge 
in CaCuOz^i have shown evidence for low energy two 
magnon scattering. Since in the near future this low en- 
ergy window will open for inelastic x-ray studies, it will 
afford an opportunity to study the dynamics of magnon 
excitations via the charge degrees of freedom which will 
complement neutron and Raman scattering studies. 

Lorenzana and Sawatzkyi 2 ^ investigated the proper- 
ties of phonon-assisted absorption in 2D Heisenberg anti- 
ferromagnets, and investigated the bi-magnon spectrum 
probed by infrared conductivity. In this paper we extend 
these calculations to present a theory of inelastic x-ray 
scattering of the two-magnon response in a Heisenberg 
antiferromagnet. In particular we show that the momen- 
tum and polarization dependence can provide detailed 
information on the nature of magnon-magnon interac- 



tions in parent insulating cuprate compounds. 

We remark at the outset that we neglect specific res- 
onant matrix elements relevant to the RIXS process in- 
cluding transitions resulting from the creation of the core 
hole, and work in the restricted model of the half-filled 
single-band Hubbard model on the square lattice to cap- 
ture scattering via magnon creation. Thus we neglect 
specific pathways of charge excitation, such as the Cu 
K-edge ls-4p transition, and how double-spin flips may 
occur near the site where the core hole is created. Such 
kinematic details are indeed important to determine ac- 
curately the RIXS intensity as well as the proper sym- 
metry and polarization of spin excitations which may be 
probed by specific x-ray transitions. However, in order 
to obtain a preliminary understanding of how the two- 
magnon response can be probed and how polarization 
may enter, we first simply focus on the evolution of the 
two-magnon response, far from any specific resonance, 
for nonzero photon momentum transfers q. 

Transitions via light scattering can be created via 
dipolc or multipole matrix elements involving states 
within the conduction band or out of the valence band. 
These transition may be selected by orienting incident 
and scattering polarization light vectors, e.; and e^, re- 
spectively. The photons entering in the scattering process 
are represented by {uJij, ^ij, &i,f) where indices i and 
/ represent the incoming and outgoing photon, respec- 
tively. Since we are interested in the insulating phase, 
the scattering of light is induce d by th e interband part of 
the operator. Following Refs. 
momentum transfer q, q = k; 
for different polarization geometries. The interband part 
of the scattering matrix clement is given by: 
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where v represents states out of the lower 
Hubbard band, and Jk is the current op- 
erator J k = Ep, CT ^T C p+|^c p _| !(T , with 
e p = —It [cos(p x a) + cos(p y a)] for a square lattice 
with lattice constant a, and nearest-neighbor hopping t. 
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The current may be expressed as: 



p<5 -ik f -(r+<5/2) 
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= E r ,i,(j 6* Jk f (r, S), 

where is the projection of the polarization vector along 
the neighbor direction 6. In the insulating state, \i) and 
/) are single occupied states, whereas \v) are excited 
states consisting of a doubly occupied state at site r + 8. 
Thus the energy of the excited state is roughly given by 
e„ ~ [/ + £«, and the intermediate states may be collapsed 
using the identity \ — S, • S,- = ^cj „c 



cl „/Cj )ff ', which 



simply means that the exchange process (and therefore 
a doubly occupied site) is allowed only if spins a and a' 

are opposite^ 4 -. Finally, we recast (f\M r \i) in Eq. 1 via magnon-magnon interactions: Io(q,cj) 



operator which, being proportional to Ji, commutes with 
7i, giving no Raman scattering in that channel. 

One also finds Oa x ,Bi vanish for q = Q = (ir, ir) and 
(2) symmetry related points for both Ai g and B\ g . This is a 
consequence of including only nearest neighbor spin in- 
teractions in the Heisenberg model 21 -. While in our case 
the response vanishes for the antiferromagnetic reciprocal 
lattice vectors, if one includes longer range interactions, 
then this restriction may be lifted 1 ^. Thus the x-ray Ra- 
man response for these wavevectors may provide a win- 
dow to sample the role of longer-range spin interactions. 

The scattering intensity is proportional to {,f\Os lg \i} 2 , 
while satisfying energy conservation u> = t^k + ^k+q- We 
first focus on the B\ g channel, where the two-magnon 
scattering is most prominent for q — 0. The left panel of 
Fig. 1 shows the B\ g bare scattering intensity, neglecting 



the following scattering operator: 

0(q) = 8t 2 EM eiq ' (r+5/2) (^e/)(^e,) 
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For crossed polarizations, e^/ = ± y), transform- 

ing as Big, and for parallel polarizations along x and y, 
transforming as A\ g , we obtain the following expression : 
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with Pai,bi(^)— 1> for 8 = an d equals 1, -1 for 5 = ay 
for Ai g and Bi g , respectively. This is a finite momentum 
generalization of the usual Loudon-Fleury light scatter- 
ing operator, and for q = the above expressions give 
the standard Raman result o 13 ! 14 . We note that in gen- 
eral subtractions of spectra for different polarization ori- 
entations arc needed in order to fully extract symmetry 
deconvolved spectr 

Henceforth, we restrict ourselves to the half-filled Hub- 
bard model with t/U small, and focus on its spin-i an- 
tiferromagnetic Heisenberg representation with Hamil- 
tonian H. We first investigate the q-dependent in- 
clastic response via the spin-wave (SW) approxima- 
tion. We proceed as usua ^ 14 ' 15 ' 16 i 17 i 18 i 19 i 20 and ex- 
press H in its SW representation, Hsw = Cst + 

E k ^k (a^Qk + PiPk) , with uj k 

magnon dispersion, and 2m 2 . — 1 = 

and 7k = (1/2) [cos(k x a) + cos(k y a)]. Performing 
the same transformations to the scattering operators 
®A lg ,B lg we obtain: 



JSZ^/l - 7 2 the 



cjk+q — w )^k q' w hcrc 6k. q is the term in curly brackets in 
Eq. 5. The Io(q — 0,oj) intensity recovers the standard 
Raman response, which for the non-interacting case has 
a peak at u — AJ due to the large magnon density of 
states at k = (ir, 0) projected out by the Bi g operator. 
For nonzero q, two magnons are created with wavevectors 
k of different magnitude and direction, leaving behind a 
reorganized spin configuration, and the response is given 
by convolving two magnon density of states at different 
k separated by q. We note that while the response is 
still governed by ~ 4J for finite q due to the magnon 
density of states being largest at k = (ir, 0), the intensity 
is highly suppressed as q approaches Q. 
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FIG. 1: Spin-wave intensity spectrum without (7o(q,w), left) 
and with (7R,pA(q, w), right) magnon-magnon interaction. 
The color intensity scale is 2.2 times smaller for the right 
panel than the left. 



OA l9 ,B lg cxE k /3R +q {- [cos(^)±cos(^)] 

X (Uk«k+q + WkMk+q) + (WkWk+q + V k V k+ q) 

x [cos(k x a) cos(3|^) ± cos(k y a) cos(^)] } , 

(5) 

where we have neglected the prefactor in Eq. 0J For A\ g 
and q = 0, we recover the familiar form of the Raman 



Magnon-magnon interactions need to be in- 
cluded to describe more realistically the local spin 
rearrangcmen t 9 i 14 ' 15 i 16 ' 17 i 18 ' 19 i 20 . In particular, interac- 
tions lead to a reduction of the peak frequency for B\ g 
Raman to 2.78J, where J, as shown by Singh et al££, 
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is also renormalized by quantum corrections. Many 
diagrams contribute to the magnon-magnon vertex 
corrections, and a subset has been investigated for q = 
Raman scattering^. In that case, magnon-magnon 
interactions lower the relevant energy scale from 4J to 
3 J due to the local breaking of six exchange bonds for 
two neighboring spin flips. For finite q considered here, 
the two magnons are created with net momentum q 
which distribute the spin flips over longer length scales. 
Since magnon-magnon vertex corrections are expected 
to weaken as the spin arrangement occurs at larger 
length scales, we approximate the renormalized response 
by a generalized random-phase approximation (RPA) 
form / RPA (q,u)) ~ q^pffe^) , where the function 
J(q) = [cos(q x a/2) + cos(q y a/2)]/2 is taken to recover 
the Raman form at q = and properties of the solution 
to the Bethe-Salpeter equation for finite q. A proper 
treatment of magnon-magnon interactions for all q in the 
SW framework is a topic of future research. As shown 
in the right panel of Fig. [TJ the magnon interactions 
for larger q bring the peak down to ui = 2.78J for B\ g 
q = Raman, but the general weakening of the magnon 
interactions at larger q gives back the unrenormalized 
response at larger q with, in particular, the response still 
vanishing at q = Q. On the other hand, the dispersion 
of the peak changes dramatically when magnon-magnon 
interactions are included, where the peak hardens at 
finite q from the r-point, q = 0. 

a R A, + B. 
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FIG. 2: 16-site cluster : exact diagonalizations for Ai g and 
Big polarizations, to be compared with spin-wave results. 
The spectra are shown for the independent momenta q al- 
lowed for the N = 16 site cluster (given on the side of right- 
most panel). The relative intensity of the curves is given by 

<o(q) 2 >- 



The results above differ slightly from those obtained 
in Ref. fl2| for infrared absorption from multi-magnons. 
While the dispersion of the peaks is similar to our results, 
Ref. [l2j found a bi-magnon response which is dominated 
by a sharp peak with a small bi-magnon lifetime for mo- 
mentum transfers (7r,0). Since the form factors are dif- 
ferent for infrared and polarized x-ray Raman measure- 
ments, the bi-magnon spectrum has different projections 



for the results presented here, favoring q = for bi- 
magnon in Big channel, and a growth of A\ g component 
with increasing q. Thus the Ai g response is more similar 
to the infrared one, as it weights out similar regions in 
the magnetic Brillouin zone (BZ). We note that phonon- 
assistcd bi-magnon response for x-ray Raman, considered 
in the same vein as Ref. fl2j . would also have similar 
polarization dependent form factors as those considered 
here, coming from the symmetry classification of the ac- 
tive phonon modes involved in the scattering. 

In order to explore the semi-quantitative validity of the 
SW/RPA results, we investigate the q-dependent inelas- 
tic x-ray scattering for the Heisenberg spin Hamiltonian 
7i, using an exact diagonalization approach, consider- 
ing both the Ai g and Bi g channels. Although exact- 
diagonalization is limited because of the prohibitive size 
of the Hilbert space, Gagliano and Balseiro^ demon- 
strated that it is a powerful technique that allows to 
compute the dynamical quantities easily. Noteworthy 
in the context of numerical investigations of Raman 
spectra, Sandvik et al^ showed that even though the 
spectra obtained by the Lanczos method are extremely 
size-dependent, they are of direct relevance to test the 
response calculated from a spin- wave analysis. As in 
Ref. [24j| , the spectra are evaluated by computing the 
continued fraction: 
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with z = ui+ie+Eo, where Eq is the ground-state energy, 
e is a damping factor. 

We first compute the ground-state energy. For N = 16 
sites we have Eo/(NJ) « —0.70178020 and wave- vector 
l^ff.s.)' We then begin evaluating the continued fraction 
with the starting state: 

l*(q)> = 



0(q)|* 9 . a .) 



v^.lon-qWq)!*^.) 



(7) 



Here, compared to Raman scattering, care has to be 
taken in computing Eq. (J6j) : the state |$(q)) being in 
a different subspace at nonzero q than the ground-state, 
the matrix elements of H. to compute the continued frac- 
tion have to be expressed in the corresponding subspace. 
The results are summarized in Figs. [5J3] 

For the case of B\ g scattering, we recover prior results 
for q = OH, while we see that the peak disperses to higher 
energies, approaching 4J both for q along the BZ diago- 
nal as well as along the axes. In light of the RPA results 
in the right panel of Fig. 1, we attribute this dispersion 
as a weakening of magnon interactions at larger q. At 
the same time, the overall intensity diminishes for larger 
q due to the form factor appearing in Eq. 

Since exact diagonalizations deal with small clusters, it 
is very difficult to make a quantitative finite-size scaling 
analysis for the spectral shaped. Therefore we repeat the 
calculation for the Bi g polarization with a 20-site clus- 
ter, and plot the results in Fig. [3] It is well known that 
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FIG. 3: 20-site cluster : exact diagonalizations in the Bi g 
channel. The inset is a schematic of the investigated q-points. 



the 16-site cluster has additional symmetries (hypercubc- 
like)^ and it therefore important to check for similar q- 
points in the 20-site cluster. Due to the different shape 
and boundary conditions of this cluster, we are not able 
to explore the same momenta in the BZ. However direct 
inspection of the q points (0,0), (tt,0) and (tt,tt) reveal 
that many of the features are qualitatively similar: the 
q=0 peak is at lower frequency than (7r,0), and the in- 
tensity weakens at finite q = |q| as in the 16-site case 
(compare with the middle panel of Fig. [2]). In addition, 
the dispersion along the path in the q-space is of the 
same order ~ J. In the RPA calculation, the minimum 
frequency of the main peak is at (0,0). For the com- 
patible q-points, the results are consistent with the RPA 
approach as well (same qualitative dispersion and mini- 
mum). We therefore believe that exact diagonalizations 
give a correct qualitative description of the response. 

For Ai g , a peak appears for finite q2i. Generally the 
peak disperses towards 4J for larger q, following essen- 
tially the Bi g spectrum. The Bi g intensity falls with 
increasing q while at the same time the Ai g spectral in- 
tensity grows. In fact, we note that the intensity for A\ g 
approaches that of the B ig spectrum for q = (tt, 7t/2), 
where the scattering operators in Eq. dl]) are the same. 

Thus, the overall intensity in an unpolarized measure- 
ment, given by Ia 1b + Ib 1s as shown in Fig. [51 is domi- 



nated by the Bi g channel for small q, and a mixture of 
Ai g and Bi g for larger q —. The main results are thus 
(1) for small q, the two-magnon Raman peak rapidly dis- 
perses upward in energy, both for q along the BZ diagonal 
and BZ axes. At the same time, a contribution from the 
Ai g channel develops, as the intensity of the Bi g falls 
with increasing q. We note, however, that the problems 
of describing the two-magnon profile for q = 0, which in- 
cludes the differences between Ai g and B\ g intensities, 
as well as the width of the lineshapes, remains an issue 
here as well (see, e.g., Ref. @). Thus more detailed treat- 
ments including further spin-spin interactions, ring ex- 
change, and possibly spin-phonon coupling, would need 
to be incorporated in low energy inelastic x-ray scattering 
as well. Besides affecting the overall lineshape, the be- 
havior at different momentum points, such as the special 
point at (tt, tt) may change. Thus, the detailed momen- 
tum dependence of the spectra may be able to provide 
important information of the types and extensions of spin 
interactions in antiferromagnets. 

Finally we remark on our results in terms of the cur- 
rent capabilities of RIXS experiments. Empirically, from 
fitting the two-magnon position and spin wave velocity, 
J in the cuprates is estimated to J « 0.13 eV. Thus the 
two-magnon Raman peak would lie generally obscured in 
the elastic line at currently available resolution, which is 
several orders of magnitude larger than the inelastic con- 
tribution. It is a possibility that the two-magnon contri- 
bution would emerge from underneath the elastic line at 
larger momentum transfers q as the peak in x-ray Raman 
disperses to 4 J. If the energy resolution can be enhanced, 
we propose that perhaps a clear two-magnon contribution 
will become visible. However, a proper treatment of the 
resonant matrix elements needs to be considered, which, 
while not affecting the dispersion, may change the rela- 
tive intensity of the spectra at specific q. This remains 
a topic for future research^ 

Note added in proof: After completion and submission 
of our work we became aware of a closely related preprint 
by Donkov and Chubukov ( cond-mat/0609002[ ) The re- 
sults are similar at small momenta but they differ for 
q = (tt, tt) due to a slight difference in the form of their 
scattering operator. 
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